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Tensor effects on the A'^ = 40 gap evolution of A'' = 40 isotones are studied by employing the 
Skyrme-Hartree-Fock-Bogoliubov (SHFB) and relativistic Hartree-Fock-Bogoliubov (RHFB) theo- 
ries. The results with and without the inclusion of the tensor component are compared with the 
experimental data. When the tensor force is included, both of the two different approaches are 
found to give the same trend and agree with the experimental one, which indicates the necessity of 
introducing the tensor force in the evolution of A'^ = 40 subshell and on the other hand the reliability 
of the methods. Furthermore, it is shown that the gap evolution is primarily determined by the 
corresponding tensor contributions from vr and p-tensor coupling in the relativistic framework. 
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With worldwide development of Radioactive Ion Beam 
(RIB) facilities, the nuclei far from the valley of /3- 
stability, i.e., the so-called exotic nuclei, have become 
more accessible in recent years [1] . The investigations on 
the shell structures of such nuclei have formed a new 
frontier in modern nuclear physics due to some novel 
phenomena discovered timely. The typical examples are 
the emergence of some new magic numbers |21 [5] and 
the quenching of some conventional shell closures [H |3] 
when approaching the isospin limits of the realistic nu- 
clei. As suggested by Otsuka and the collaborators, the 
tensor force may be one of the crucial physical mecha- 
nisms in modeling the shell evolution and the occurrence 
of magic numbers in the exotic regions [SHZI- However, 
in both non-relativistic and relativistic self-consistent 
mean-field calculations, the relevant tensor components 
were dropped up to very recently [51113). 

Since the magic character of ^®Ni was proposed in the 
early eighties by Bernas et al. jM], the shell closure of 
iV = 40 has been devoted more and more efforts from 
both experimental and theoretical sides due to the sig- 
nificant roles played by the nuclei of N/Z = 40 ~ 60 
involved in the r- and rp-process paths [TS]. Although 
the existence of the iV = 40 shell gap is well recognized 
in ^^Ni, away from that the uncertainties have remained 
for many years on the persistence of the shell closures. 
Until recently several experiments of /3-decay indicate a 
rapidly weakening magicity of iV = 40 away from ^^Ni 
|15H18j . whereas such trend cannot be properly repro- 
duced by shell model and mean-field calculations when 
the tensor components of the nuclear interaction are ex- 
cluded dSfrai. 

In this Brief Report, the isospin evolutions of sub- 
shell = 40, extracted from the self-consistent calcu- 
lations of A^ = 39,40,41 isotones, are studied within 
the non-relativistic and relativistic mean field theo- 
ries, specifically the non-relativistic Skyrme-Hartree- 
Fock-Bogoliubov (SHFB) [20], the relativistic Hartree- 



Bogoliubov (RHB) [2r,'22l, and relativistic Hartree-Fock- 
Bogoliubov (RHFB) j23j theories. Particular attention 
is payed to the effects of tensor couplings embedded in 
the neutron-proton interaction. In both non-relativistic 
and relativistic cases, the pairing correlations are treated 
within the Bogoliubov scheme such that the continuum 
effects can be naturally involved, which is meaningful to 
provide self-consistent calculations for the selected iso- 
tones. 

Within the SHFB theory, the tensor interaction is in- 
troduced as [50] 
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where k = i — ^2)/(2«) acts on the right and k' = 
2)/(2i) acts on the left. The coupling con- 
stants T and U denote the strength of the triplet-even 
and triplet-odd tensor components, respectively. In fact 
it has been demonstrated that the tensor force [see eq. 
([1])] plays an important role in determining the evolution 
of the nuclear shell structure [231 US] . 

Notice that the tensor force of Eq. ([!]) is in a non- 
relativistic form, inconsistent with the Lorentz covariant 
principle. Within the relativistic framework, such tensor 
correlations in T = channels can be naturally intro- 
duced with the inclusion of the tt and p-tensor couplings, 
e.g., in the relativistic Hartree-Fock (RHF) mean field 
[T21[55]. Within RHFB, the tt and p-tensor couplings are 
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adopted as [12] 

drdr'4,{r)^j{r')T^{r,r')D^{r,r')^j{r')4,{r), 

(2) 

where (j) denote the pseudo-vector tt and tensor p cou- 
phngs, D^{r,r') are the Yukawa propagators, and the 
interaction matrices T^{r,r') read as, 

r.(r,r') = - J^(/,r757M9^)r • iUrj^J.dnr' , (3a) 

rj(r,r') ^^(f^a,,rd% ■ {fpO^'fdi)r'. (3b) 

In the above equations M(to^) denotes the mass of the 
nucleon (tt meson) and fjrifp) is the coupUng constant 
of 7r(/9) meson which is exponentially density-dependent 
[T2] . Here what should be clarified is that the tensors ex- 
pressed above are associated with the Lorentz rotation, 
different from the quantities defined by the space rota- 
tion, e.g., the ones in Eq. ([I]). In fact it is also illustrated 
that the tt and p-tensor couplings, particularly the ten- 
sor components reduced with non-relativistic limit, play 
a crucial role in the self-consistent description of shell 
structure and corresponding evolution [T2 HSl US] • 

To clarify the tensor effects in modeling the shell evo- 
lution of = 40, different effective interactions are em- 
ployed in the present work. For the SHFB theory, the 
effective interaction SLy5 [27] and the one with tenor 
components added perturbatively (marked by SLy5-t-T) 
are adopted [55]. Besides SLy5, another Skyrme force 
T3I and the one with the tensor force (T3I-1-T) are 
also utilized for comparison. In the RHFB calculations 
we choose the effective interactions PK02, PK03 and 
PKAl, comparing to those by RHB with DD-ME2 [55] 
in which the Fock terms are dropped. It should be men- 
tioned that the effective Lagrangian PK02 does not con- 
tain TT or tensor p couplings, whereas the former is in- 
cluded in PK03, and both are taken into account by 
PKAI. 

Aiming at the isospin evolution behavior of sub-shell 
iV = 40 close to the neutron drip line, we select six 
neutron-rich isotones from Z = 20 to 30, which are spher- 
ical or nearly spherical as predicted by both macroscopic- 
microscopic model [30] and self-consistent mean-field the- 
ory [19]. Experimentally it is hard to measure the sub- 
shell TV = 40 directly, especially for the exotic isotones, 
e.g., ^°Ca. As the experimental references, such shell 
gaps can be approximated by the difference of single- 
neutron separation energies Sn around the magic or semi- 
magic neutron number [30l |3T] 

i^gap(^, Mn) ^SniZ, N^) - S„{Z, + 1) (4) 

In fact our RHFB calculations also show that the evalu- 
ations of sub-shell iV = 40 by Eq. ^ agree with the en- 
ergy gaps between the corresponding single-particle lev- 
els to certain precisions (~200 keV). For convenience we 
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FIG. 1: (Color online) The neutron gap at = 40 of A/' = 40 
isotones as functions of the proton number Z. The results 
are calculated by the non-relativistic (a) and relativistic (b) 
mean-field theory without and with the tensor component. 
Experimental values approximated by Eq. Q are also shown 
for comparison. See the text for details. 



choose the later as the theoretical sub-shells of = 40 
in the following discussions. 

In Fig. [l]are shown the sub-shell gaps of = 40 along 
the isotonic chain of = 40, extracted from the cal- 
culations of the non-relativistic [plot (a)] and relativistic 
[plot (b)] mean field theories. As the reference, the exper- 
imental data [32] extracted by Eq. Q are denoted by the 
black balls. It is found that the experimental Eg^p is con- 
sistent with the /3-decay studies by Mulleller et al. [TT] 
that the semi-magicity of A^ = 40 becomes weak when 
moving away from <^*^Ni. From Fig. [ija) one can clearly 
see that the experimental trend cannot be reproduced 
by the Skyrme forces SLy5 and T31 without tensor com- 
ponents. In contrast, when the tensor force participates 
in the neutron-proton interaction, identical systematics 
with the data can be provided by the calculations with 
SLy5-|-T and T31-I-T. Such distinct improvement can be 
well understood from the nature of tensor force '5]. In 
the SHFB calculations, the sub-shell of A^ = 40 is de- 
termined by the energy gap between the neutron orbits 
(or i/pi/2) and 1^59/2, i.e., the ones with j< =^-1/2 
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FIG. 2; (Color online) Detailed contributions of the neutron 
gap at N — 40 along A'^ = 40 isotones from the isovector 
TT and tensor p couplings in comparison with those from the 
other channels associated with the kinetic energy, cr-scalar, cj- 
vector, p-vector couplings, and the rearrangement term. The 
results are extracted from the RHFB calculations with PKAl. 



and jy =1 + 1/2, respectively. Along the isotonic chain 
from ^°Ca to ^^Ni, the valence protons will gradually oc- 
cupy the proton orbit 7rl/7/2 with = I + 1/2, which 
presents repulsive tensor couplings with the orbit h'gg/2 
(j> =1 + 1/2) and attractive ones with vf5/2 (or vpi/2)- 
As a result the semi-magicity of iV = 40 are remarkably 
weakened by the tensor force when moving away from 
^^Ni, in accordance with the experimental reference data 
[32] as well as the /3-decay studies [17]. As compared to 
the SHFB calculations excluding tensor components it is 
then well demonstrated for the crucial role played by the 
tensor force in modeling the shell evolution [25, 28 . 

Similar systematical agreements and discrepancies are 
also found between the relativistic calculations and the 
experimental data. As seen from Fig. [l]jb) the calcu- 
lations of RHB with DD-ME2 and RHFB with PK02 
present opposite isospin evolution behavior to the ex- 
perimental trend, in which the isovector tensor compo- 
nent TT-coupling is not taken into account. In fact due 
to the lack of isovector tensor couplings similar incon- 
sistences are also found in other RHB calculations, e.g., 
with PKDD [33 . Compared to those excluding isovector 
tensor components, the RHFB calculations with PKAl 
present identical trend as the data, which is mainly due 
to the effects of tensor components embedded in isovec- 
tor TT and tensor p couplings, particularly the former one. 
As shown in Fig. [2] only the contributions from tt and 
tensor p-couplings are consistent with the experimental 
trend while the other channels present fairly remarkable 
cancellation with the referred behaviors. Consequently 
the calculations with PKAl lead to fairly weak isospin 
dependence on the shell evolution when comparing with 
the non-relativistic results and the experimental trend. 
In fact the tt contributions extracted from the calcula- 



tions with PK03 also present consistent trend with the 
experimental data (see Fig. [3| while too strong cancel- 
lations are contributed by the other channels, leading to 
opposite isospin evolution in total [see Fig. [T](b)]. 

It is interesting that the non-relativistic SHFB calcu- 
lations with the tensor force and the relativistic RHFB 
ones including the tt- and p-tensor couplings favor the 
same trend. Different models may validate each other, 
which to a large extent indicates the reliability of the the- 
oretical models and thus one can obtain more compelling 
results. In addition, recent experimental |34[ I35| and 
theoretical |36j studies seem to indicate that a quantum 
phase transition from the spherical to deformed shapes 
takes place near = 40 for Cr and Fe isotopes, referred 
to as a "new island of inversion" . Similar to the one of 
N = 20 [37HS2, the island of inversion of = 40 could 
be explained by the weakening semi-magicity of A^ = 40 
in the neutron-rich region, which can only be described 
appropriately with the inclusion of the tensor force as 
one can see from Fig. [l] and [2] 
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FIG. 3: (Color online) Detailed contributions from the tensor 
and central parts of pion (upper panel) and p-tensor (lower 
panel) couplings to the energy difference AE = E^ig^^^ — 
between the neutron canonical single-particle states 
i'lgg/2 and 1^1/5/2 along the isotonic chain ol N = 40. The 
results are extracted from the RHFB calculations with PKAl 
(filled symbols) and PK03 (open symbols). 

To further clarify the contributions from the tensor 
coupling in RHFB theory, we shown in Fig. [3] the contri- 
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butions from the tensor and central parts of pion (upper 
panel) and p-tensor (lower panel) couplings to the evolu- 
tions of the energy difference ^Egj = E^ig^^,^ — E^if^^^, 
which corresponds to the = 40 shell gap in PKAl cal- 
culations along the isotonic chain. Here the tensor part of 
/9-tcnsor coupling is obtained from similar non-relativistic 
reduction as the tt pseudo- vector coupling (SHI SO] ■ For 
the isovector tt coupling which in total presents consis- 
tent results with the data (see Fig. [T]), its tensor part 
dominates the isospin dependence whereas the central 
part presents nearly constant isospin evolution in both 
PKAl and PK03 calculations. On the contrary, both 
the tensor and central parts of p-tensor coupling provide 
almost equivalent contributions to the isospin evolution 
of ^Egf, both are opposite with the tt effects. It is also 
due to such remarkable cancellation between the tt and 
/9-tensor couplings that PKAl provides consistent while 
fairly weak isospin-dependent behavior as referred to the 
data. 

In summary, the tensor effects in modeling the semi- 
magicity of = 40 are studied within the non-relativistic 



Skyrme-Hartree-Fock-Bogoliubov (SHFB) theory and 
the relativistic Hartree-Fock-Bogoliubov (RHFB) theory. 
The calculations with and without isovector tensor com- 
ponent compared with the experimental data indicate 
that the tensor forces play a crucial role in reproducing 
the isospin evolution of the = 40 subshell. It is also 
interesting that, when the tensor force is included, both 
present non-relativistic and relativistic methods yield the 
same trend for this shell gap evaluation. These two dif- 
ferent models validate each other and lead to a more con- 
vincing results. In the RHFB calculations, both 7r-tensor 
and /9-tensor couplings present significant contributions 
in the T = NN interaction, which eventually determine 
the gap evolution of the iV = 40 subshell, in agreement 
with the experimental trend. 
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